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Abstract 

Using the ETFSI (extended Thomas-Fermi plus Strutinsky inte- 
gral) method, we have calculated the fission barriers of nearly 2000 
exotic nuclei, including all the neutron-rich nuclei up to A = 318 that 
are expected to be relevant to the r-process, and all superheavy nuclei 
in the vicinity of N = 184, with Z < 120. Our calculations were per- 
formed with the Skyrme force SkSC4, which was determined in the 
ETFSI- 1 mass fit. For proton-deficient nuclei in the region of N = 
184 we find the barriers to be much higher than previously believed, 
which suggests that the r-process path might continue to mass num- 
bers well beyond 300. For the superheavy nuclei we typically find 
barrier heights of 6-7 MeV. 
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1 Introduction 



The r-process of stellar nucleosynthesis depends crucially on the masses and 
fission barriers (among other quantities) of nuclei that are so neutron-rich 
that there is no hope of being able to measure them in the laboratory (see 
Refs. [0, U for reviews discussing the nuclear data required for an under- 
standing of the r-process). It is thus of the greatest importance to be able to 
make reliable extrapolations of these quantities away from the known region, 
relatively close to the stability line, out towards the neutron-drip line. Un- 
til recently the masses and barriers used in all studies of the r-process were 
calculated on the basis of one form or another of the liquid- drop (let) model 
(LDM). However, in an attempt to put the extrapolations on as rigorous a 
footing as possible we have developed a mass formula that is based entirely 
on microscopic forces, the ETFSI-1 mass formula || [|, |J ||, [7|. Calculations 
of the r-process using the ETFSI-1 masses have already been performed || § . 
but they are incomplete in that fission had to be neglected, barriers not yet 
having been calculated in the ETFSI model. Of course, recourse could have 



been made to the extensive barrier calculations |K| [U]] based on the LDM 
that had been used in earlier r-process studies, but it would have been incon- 
sistent to use one model for the masses (required for the neutron-separation 
energies S n and the beta-decay energies Qp), and another for the barriers. 

Here we remedy this deficiency by presenting the results of ETFSI-method 
calculations of the fission barriers of all of the nearly 2000 nuclei lying in the 
region of the (N,Z) plane shown in Fig. 1. This region covers the range 
84 < Z < 120, and extends from moderately neutron-deficient to extremely 
neutron-rich nuclei, including thereby not only a large fraction of the nuclei 
whose barriers are known experimentally, but also all nuclei with A < 318 
whose barriers can reasonably be expected to be relevant to the r-process. 
Towards the upper limit of our range of Z values the r-process path lies at 
much higher values of A (if it has not already been terminated by fission), 
and we restrict ourselves to nuclei lying close to the stability line in the long- 
expected "island" of stability that is now becoming experimentally accessible 
0,11,11. 

The ETFSI method is a high-speed approximation to the Skyrme-Hartree- 
Fock (SHF) method, with pairing correlations generated by a 5-function 
force, treated in the usual BCS approach (with blocking). There are two 
parts to the total energy calculated by the ETFSI method, the first con- 
sisting of a purely semi-classical approximation to the SHF method, the 



2 



full fourth-order extended Thomas-Fermi (ETF) method, while the second 
part, which is based on what we call the Strutinsky-integral (SI) form of 
the Strutinsky theorem, constitutes an attempt to improve this approxima- 
tion perturbatively, and in particular to restore the shell corrections that 
are missing from the ETF part. The way in which we extended the ETFSI 
method, developed originally for the calculation of ground-state energies, to 



the large-scale calculation of fission barriers is described in Ref. JT^|, which 
should be consulted for all details of our calculational methods. However, 
it is important to recall here that our fission paths are optimized with re- 
spect to the the elongation parameter c, the necking parameter h, and the 
asymmetry parameter at (see Ref. |15l , and in particular its App. A, for the 



definition of these quantities, and a description of the nuclear shapes allowed 
by our parametrization in the case of fission; note that as in that reference 
we are assuming axial symmetry). 



In Ref. |15j we calculated all the barriers that have been measured in 



nuclei with Z > 81. The comparison of our results with the available data 
showed that if the primary, i,e., highest, barrier is lower than 10 MeV the er- 
ror never exceeds 1.4 MeV, with either sign being possible, while for primary 
barrier heights lying between 10 and 15 MeV we probably overestimate the 
height by about 1.5 MeV; barrier heights in excess of 15 MeV are overesti- 
mated by 2 MeV or more. Since the primary barrier is the only one that is 
really relevant to the r-process we see that our method is sufficiently accu- 
rate for the calculation of barriers that are low enough to be of interest to 
the r-process (see Section 3 below). The calculations here are performed in 
exactly the same way, and in particular we use the same force parameters, 
set SkSC4, which, it should be noted, were determined entirely by the mass 
fit ||, and have in no way been modified for the barrier calculations. 



2 Results 

We summarize the results of our calculations in Fig. 1, and in Table [I], 
where for each nucleus we present the primary barrier. In order to keep its 
size within reasonable bounds, Table [l| does not show the three deformation 
parameters of the barriers, c, h, and a. However, these parameters can be 
found on our web site (jlittp:/ /www-astro. ulb.ac.be/| ), which gives also the 
value of the secondary (i.e., second-highest) barrier, when present, as well 



as its deformation parameters. As mentioned in Ref. [15], it is possible to 
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Figure 1: Schematic representation in the (N,Z) plane of the primary fission 
barriers displayed in Table [I]. The energy bins are described in the figure. 
Thin solid lines represent r-process paths for constant values (in MeV) of S a 
(defined in text). The neutron-drip line is shown by the thick solid line 
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distinguish without ambiguity between "inner" and "outer" barriers from 
their value of c, each of the two highest barriers belonging to a different 
category. A typical situation is illustrated in Fig. 2 for the uranium isotope 
chain. For smaller values of Z the separation between primary and secondary 
barriers, expressed in terms of c, tends to increase, and the external barrier 
becomes the primary one for all isotopes, as it is already in Fig. 2 for large 
neutron excesses (this can be explained by the fact that nuclei with smaller 
fissility need higher deformation to fission). For larger Z this separation 
generally decreases and the inner barrier is always the highest, while the 
outer one progressively disappears, so that for Z > 100 only one barrier (the 
inner one) is left. We label the primary barriers displayed in Table [I] with 
the superscript % or o to show whether they belong to the inner or to the 
outer category, respectively. 

The general trends of the results of Table [IJ are displayed in Fig. 1 as 
well as by the curves labelled "ETFSI-1" in Figs. 3-5, where we show three 
isotope chains, Z = 84, 92, and 100, respectively (note that in Fig. 3 the 
ETFSI-1 barriers for A = 207-212, not included in Table [1|, are taken from 



Ref. 15 



The curve labelled "ETF" denotes the results we find for pure ETF cal- 
culations of the barriers without either shell or pairing corrections. In these 
latter calculations we search anew the positions of the ground states and 
saddle points, as defined by their three deformation parameters, and find 
different deformations for these points than obtained in the ETFSI calcu- 
lation. We emphasize that if shell and pairing corrections had simply been 
subtracted from the ETFSI saddle-point and ground-state energies, the ETF 
curve would not be at all as smooth as it is in Figs. 3-5. 

Also seen on these graphs are the results of the LDM calculations of 
Howard and Moller (labelled "HM") [PT(] , and the results of Myers and 
Swiatecki (labelled "MS" ) |16[ , based on zeroth-order Thomas- Fermi calcula- 
tions Both these latter sets of results include shell corrections calculated, 



in one way or another, by the Strutinsky method. Actually, Ref. [16j gives 
just a smooth formula representing the main trends of Thomas- Fermi barrier 
calculations, which we reproduce here as the curve "MS.0". We have con- 



structed the curve "MS" ourselves, following the prescription of Ref. [TT>| . 
i.e., we have added the shell corrections quoted in that paper (see their Ref. 
[4]) to the ground-state energy only, assuming there are no shell corrections 
at the saddle points (this is their "topographic" theorem). 
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Figure 2: (a) Inner (solid curve) and outer (dot-dashed curve) barriers for 
the Z = 92 isotopic chain. The (single) ETF barrier is also shown (dashed 
curve), (b) The corresponding elongation parameter c. 
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Table 1: Heights (in MeV) of primary barriers ( l denotes inner, outer - see 
text) 
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3 Discussion 



3.1 Comparison with experiment 



The solid circles in Figs. 3-5 denote the few experimental barrier heights 
that have been measured for these isotope chains [HJ. In general, all three of 
the calculations "ETFSI-1", "HM", and "MS" are in good agreement with 
the data, although in the case of Z = 84 the "ETFSI-1" barriers are between 
2 and 4 MeV too high. This discrepancy reflects the general tendency noted 



above and in Ref. |15|]: we overestimate all barriers higher than 15 MeV by 
2 MeV or more. Possible reasons for this are discussed in Ref. TBI; another 



point not mentioned there is that barriers will be much more sensitive than 
masses to the droplet-model curvature coefficient a cv (HJ corresponding to 
the Skyrme force. Thus it is conceivable that despite the good mass fit our a cv 
is sufficiently overestimated to lead to serious errors at extreme deformations. 
(The Thomas-Fermi calculations of Ref. |T7[ encountered the same problem, 



but resolved it by adjusting the so-called "congruence energy"; we stress that 
we do not make use of this feature in our calculations, and indeed if we had 
done so then the good agreement with experiment that we find for the lower 
barriers would have been destroyed.) In any case, these large errors are of no 
practical consequence for the r-process, since nuclei with such high barriers 
will be effectively stable against fission in this context. On the other hand, 
for all barriers lower than 15 MeV the ETFSI-1 results never disagree with 
experiment by more than 1.7 MeV, and usually by much less: for Z > 86 
the rms error is 717 KeV, while for Z > 88 it is as small as 698 KeV (this 
includes a few measured nuclei that are given in Table |2] of Ref. [pl| but not 
here) . 



3.2 Highly neutron-rich nuclei 

While all three methods agree reasonably well over the narrow experimentally 
known region, they give widely different extrapolations to the highly neutron- 
rich region. In particular, the ETFSI-1 calculations predict exceptionally 
high barriers in the region of iV = 184 for the lower values of Z, although 
as Z increases, i.e., as the stability line is approached, these barriers become 
lower, and the various sets of calculations tend to converge. Actually, a 
similar trend in the vicinity of N = 184 can be observed in the case of the 
MS calculations, although on a much lesser scale, while the HM calculations 
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Figure 3: Primary fission barriers for the Z = 84 isotopic chain, calculated 
with the various models considered in text. Neutron closed shells are indi- 
cated by vertical dashed lines, experimental values by solid circles. 
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predict no enhancement of barrier heights at all close to iV = 184. Because of 
our tendency to overestimate barriers that are high anyway, the trend in the 
ETFSI-1 results shown in Figs. 3-5 may be somewhat exaggerated. However, 
we can estimate these errors by referring to Table [I] of Ref. [|15[], and we find 
that our calculated barrier for nucleus Z = 84, N = 184 is probably no more 
than 12 MeV too high, while that of Z = 92, JV = 184 is only 2 MeV too 
high at the most. Thus the tendency we have reported for ETFSI-1 is at 
least qualitatively real, and it is certainly stronger than in the MS case. 

The rapid fall-off beyond N = 184 for all three isotope chains is indicative 
of a shell effect. In fact, the behaviour of the ETFSI-1 barriers in the vicinity 
of N = 184 must be at least partially related to the fact that the force SkSC4 
leads to a strong magic gap at N = 184 for proton-deficient nuclei, the gap 
becoming much smaller as the stability line is approached with increasing Z. 
This point is illustrated in Fig. 6, where we show the variation of the gap 
A = S 2n (Z,N = 186) - S 2n (Z,N = 184) with Z. (It must be recalled that 
we have calculated pairing in the BCS approximation, and that with the 
more realistic Bogolyubov treatment much weaker neutron shell gaps will 



be found in general for large neutron excesses [|20 |. However, unpublished 
results of Dobaczewski, discussed by Pearson et al. plf , show that this 
"Bogolyubov quenching" of shell gaps is much less pronounced for the N = 
184 magic number, and we therefore neglect it here.) We also show in Fig. 6 
the variation of this same gap for the FRDM ( "finite-range droplet model" ) 
mass formula [Ell. Since the shell corrections of the MS calculations consist 



entirely of those of the FRDM (applied exclusively to the ground state) we 
can understand why the MS barriers peak much less strongly near N = 184. 
Sensitivity of barrier heights to symmetry coefficient 

However, the very high barriers found for ETFSI-1 near iV = 184 cannot 
be the result entirely of shell effects, since the pure ETF calculations dis- 
played in Figs. 3-5 show that the macroscopic part plays a role also, there 
being a steady increase of barrier height with N, for constant Z. Referring 
to Table 0, one sees in fact that the barrier heights of these highly neutron- 
rich nuclei are strongly anti-correlated with the symmetry coefficient J of 
the model in question, or equivalently, strongly correlated with the surface- 
symmetry coefficient a ss , defined by 

a ss = (2L/K v )a sf - 9J 2 /4Q , (1) 

where L is the density-symmetry coefficient, K v the incompressibility, a s f the 



surface coefficient, and Q the surface-symmetry stiffness coefficient [19]. The 
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Figure 4: Same as Fig. 3 for Z = 92. Also shown on this figure are some 
results for forces SkSCIO and SkSC15. 
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Figure 5: Same as Fig. 3 for Z = 100 
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anti-correlation between J and a ss holds for all entries in Table and is, in 
fact, a quite general property of all models that have been fitted to nuclear 
masses, whether they are of the droplet type or are based on microscopic 
forces |H |4 . 



This result can be easily understood in terms of the following gross ap- 
proximation to the droplet model, which holds best close to the stability 
line, 

e = a v + a sf A-V 3 + ( J + a ss A-^)I 2 + a coul Z 2 A~^ 3 + . . . , (2) 

where e is the energy per nucleon and / = (N — Z)/A. We see now that for 
ground-state masses of nuclei relatively close to the stability line an increase 
in J can be roughly compensated by a decrease in a ss over a large range of 
values of A. Now the term in a ss A -1 / 3 is really a surface term, so we can 
write the fissility parameter as 

X= 2^(ljA ' {6) 

where 

a sf (I) = a sf + a ss I 2 . (4) 

Thus, while a decrease in J will be compensated by an increase in a ss as far 
as ground-state masses are concerned (at least for nuclei relatively close to 
the stability line), the result will be an increased barrier height for nuclei of 
large neutron excess I. 

This interpretation of our very high barriers as being at least in part 
related to macroscopic symmetry properties is strengthened by recalculating 
the barrier of Uig4 with force SkSCIO, a force that has been fitted in exactly 
the same way as force SkSC4, except that it is constrained to J = 32 MeV 
||25|| , rather than 27 MeV: it will be seen from Fig. 4 that this barrier is 
lowered by more than 6 MeV. At the same time, for a nucleus much closer to 
the stability line, Uuc the barrier is lowered by only 1 MeV on replacing force 
SkSC4 by SkSCIO, as we would expect for a symmetry-related effect, and we 
see that the agreement with the experimental barriers that we found in Ref. 
[IS | is more or less J-independent. However, our new barrier for Uig4 is still 



considerably higher than that given by the MS calculation, even though the 
latter has almost the same value for J (32.65 MeV). This suggests that shell 
effects are still playing a major role in our high barriers, an indication that 
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Figure 6: The N = 184 pairing gap as a function of Z for SkSC4, SkSCIO 
and FRDM. 
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is confirmed in Fig. 6, where it will be seen that the shell effects for our two 
Skyrme forces are more or less the same. 

Nevertheless, it is clear that there is a strong macroscopic symmetry 
effect, at least within the ETFSI framework, and it is obviously essential 
that we tie down much better the value of J (and thus of a ss ). Now the rms 
error of the mass fit for force SkSCIO is much worse than for the ETFSI- 1 
force, SkSC4 (a = 0.893 MeV rather than 0.736 MeV), and we now find that 
the highest value of J for which acceptable mass fits can be found in the 
ETFSI framework is 28 MeV. On the other hand, this is the lowest value 
for which the known stability of neutron matter can be assured, so we now 



regard 28 MeV as our best value for J [26]. Recent HF calculations [27] reach 
an identical conclusion, and it now seems that 28 MeV is a quite robust value 
within the very general framework of Skyrme forces. This conclusion stands 
in contradiction with the value of J 



32.73 MeV given by the FRDM 
the best droplet-model fit to masses; the almost identical value adopted in 
the zeroth-order Thomas- Fermi calculations of MS is presumably dictated by 
the requirements of self-consistency, the shell corrections of these calculations 
coming from the FRDM. Ref. p6] suggests one way in which the FRDM 
could lead to a spuriously large value of J, but in any case we have no option 
in an ETFSI-model calculation of barriers but to choose lower values of J if 
we wish to simultaneously fit masses. 



Table 2: Volume- and surface-symmetry coefficients of the forces and models 
used in this paper (see text) 



SkSC4 



MS 



HM 



SkSCIO 



SkSC15 



J (MeV) 
L (MeV) 
K v (MeV) 
a sf (MeV) 
Q (MeV) 
a ss (MeV) 
Ref. 



27.0 
-9.29 
234.7 
17.7 
75* 
-23.3 

m 



32.65 
49.9 
234 
18.63 
35.4 
-59.9 





36.5 
100 
240 
20.76 
17.0 
-159 



32.0 
55.82 
235.8 
18.11 

34 
-59.2 

wn 



28.0 
6.73 
234.9 
17.78 

56 
-30.5 



*This newly calculated value replaces the one given in Ref. 



The question now arises as to how much our barriers would be lowered 
if we changed J from the value of 27 MeV corresponding to the force SkSC4 
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with which they were calculated (and which leads to a unphysical collapse of 
neutron matter) to our newly preferred value of 28 MeV. We have accordingly 
constructed a new force, SkSC15, that has been fitted to the same mass data 
as SkSC4, but under the constraint of J = 28 MeV, rather than 27 MeV 
26| . Repeating the calculation of Ui84, we see from Fig. 4 that the barrier is 



lowered by 3.2 MeV, which still leaves it much higher than the MS and HM 
values. Moreover, this is an extreme case, and closer to the stability line the 
effect of changing J by 1 MeV will be negligible. 

The r-process 

The fact that the barriers we find for very neutron-rich nuclei are much 
higher than hitherto believed has significant implications for the r-process, 
some of which have been presented in Ref. ||28||. Here we confine ourselves to 
a few general remarks. 

For fission to occur on a timescale that is short compared to the beta- 
decay lifetimes of the nuclei found on the r-process path, and the r-process 
path thereby terminated, the nucleus in question will have to be excited close 
to or over the top of the barrier. Such excitation of a nucleus during the r- 
process can occur as the result of either neutron capture (neutron-induced 
fission) or beta-decay (beta-delayed fission), so for rapid fission to occur 
the S n or the Qp (of the parent nucleus), respectively, must be close to or 
higher than the height of the primary barrier. Now Fig. 1 displays r-process 
paths corresponding to different values of S a , defined as half the two-neutron 
separation energy S^n, as well as the neutron-drip line (all calculated with 
the ETFSI-1 mass formula). With typical values of S a lying between 1 and 
2 MeV, we see from this figure that neutron-induced fission will certainly 
not occur in the r-process below A = 318, and probably not below much 
higher values. On the other hand, referring to the ETFSI-1 mass tables 
for the relevant Qp's, we find that beta-delayed fission will certainly be 
possible for some r-process nuclei with A < 300, but the extent to which this 
occurs depends on the beta-decay strength function of the precursor being 
sufficiently concentrated towards the upper end of the spectrum. But some 
fraction of r-process nuclei will always escape beta-delayed fission, and with 
neutron-induced fission no longer being operative we now have to entertain 
the possibility of the r-process path extending to values of A considerably in 
excess of 300. This conclusion is not vitiated by a possible overestimation 
of our barrier heights, the limits of which are discussed in the previous sub- 
section. 
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3.3 Superheavy nuclei 

We have already pointed out how our barriers in the vicinity of N — 184, 
while very high on the r-process path, become lower as Z increases, but even 
as the stability line is approached they become higher again, confirming all 
previous expectations of an "island" of stability in this region. (Actually, 
isofar as one speaks of the "valley of stability" , it would be more appropriate 
to speak of a "basin of stability" than of an "island" .) Thus for nuclides in the 
range 112 < Z < 120 and 177 < N < 186 our calculated barriers are always 
at least 5 MeV, and occasionally nearly 9 MeV, high, assuring thereby a 
large measure of stability, at least with respect to fission. As for the recently 
discovered superheavies, 289 114 []TJJ and 293 118 our calculated barrier 



heights are 6.7 and 7.4 MeV, respectively. The barriers of these nuclei have 
not yet been measured, but the very fact that they are stable enough to have 
been observed at all constitutes a qualitative verification of the high barriers 
that our calculations predict. 

As for the extent to which we agree with other microscopic calculations 
in this region we note that in a SHF calculation on 288 112 Cwiok et al. find 
a barrier height of about 6.5 MeV |p9[| , as compared with our value of 4.8 



MeV. This difference is not very large, but their barrier is a double one, 
with both peaks having about the same height, whereas we have a single, 
inner, barrier. All in all, it looks very much as though this nucleus is more 
deformable in our calculation than in that of Ref . [53 , despite the fact that 



the force they use, SLy7 f30 |, has a lower surface coefficient a s f than our own 



SkSC4 (17.0 MeV rather than 17.7 MeV). The reason could be that the force 
SLy7 |30| has a low effective mass M*, 0.69M, which will lead to too low a 
density of single-particle states, and thus an unrealistically high resistance 
to deformation. (The fact that Ref. |2£| imposes reflection symmetry while 
we do not cannot be a factor, since in our calculations this nucleus turns out 
to be reflection-symmetric anyway.) 

The only other microscopic calculations in this region of which we are 
aware are those of Berger et al. who use the HF-Bogolyubov (HFB) 



method with the finite-range Gogny force. They find a barrier height of 
10 MeV for 294 112, of 11 MeV for 298 114, and of 7 MeV for 306 118; these 
results are to be compared with our own values of 6.9, 7.5, and 5.6 MeV, 
respectively. The consistently higher barriers obtained by Ref. |JT|] could be 
accounted for in the same way that we have suggested for the case of Ref. 
too low an effective mass for the Gogny force. However, an additional 
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contribution could now come from the surface coefficient, since this is higher 
for the Gogny force (20.1 MeV [[32)) than for SkSC4. (The fact that Ref. 



3lf imposes reflection symmetry while we do not can be a factor only in the 
case of 306 118, since we find the other two nuclei to be reflection-symmetric 
anyway.) 

As for the MS procedure of Ref. 



16], it cannot be applied to nuclei 



with Z > 112, since this corresponds to a limiting fissility beyond which no 
prescription is given in that paper for calculating the macroscopic barriers. 



Of course, the original TF method [|17] on which Ref. |L6| was based is just 
as applicable in the superheavy region as elsewhere, but no such calculations 
seem to have been published. 



4 Conclusions 

We have extended the ETFSI-1 mass formula, based on the Skyrme force 
SkSC4, to the calculation of the fission barriers of all nuclei that can be 
expected to play a role in the r-process of nucleosynthesis; we recall that 
the force SkSC4 gives an excellent fit to the mass data, and to the known 
primary-barrier heights that are lower than 15 MeV. The results that we 
present here are radically different from the only other such calculations that 
have been made |10|, |j~6| , in that we obtain much higher barriers for proton- 
deficient nuclei in the region of N = 184, a consequence of which is that the 
r-process path might continue to mass numbers considerably in excess of 300 
before being brought to a halt by neutron-induced fission. In view of the 
importance of this result we summarize here the reasons why we believe that 
our calculations are essentially correct. Our high barriers on the r-path are 
related both to the shell effects associated with our Skyrme force, and to our 
much lower value of the symmetry coefficient J (27 MeV as opposed to 36.5 



MeV for HM |T(J and 32.65 MeV for MS |16[]); we now examine each of these 
points in turn. 

Our shell effects appear to be fairly robust within the ETFSI framework, 
being difficult to change significantly while maintaining the fit to the mass 
data. At the same time, we recall that in extrapolating far from the known 
region of the nuclear chart out to the highly neutron-rich region that is 
relevant to the r-process, the isospin dependence of the ETFSI spin-orbit 
field conforms well to the predictions of relativistic mean-field theory |33|, |34 
adding thereby to our confidence in the reliability of the extrapolation. 
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As for the value of the symmetry coefficient J, we have found, after the 
completion of these calculations, that 28 MeV would have been a better value 
than 27 MeV, essentially because it allows us to avoid an unphysical collapse 
of neutron matter, while maintaining a high-quality mass fit. A complete 
new mass fit, ETFSI-2, is currently being undertaken with the constraint 
J = 28 MeV, and in principle when this is completed we should repeat the 
calculations of this paper with the new force determined by this fit. However, 
we have shown here that the maximum effect of the change of force will be 
to lower the barriers of the most neutron-rich nuclei by about 3 MeV, and 
that for nuclei closer to the stability line the effect will be correspondingly 
smaller. In any case, one can be reasonably sure that the ETFSI-2 fission 
barriers would lie considerably closer to the ETFSI-1 barriers presented here 



than to the barriers of either MS |16] or HM [10]. 



We have also extended our barrier calculations of nuclei with N in the 
vicinity of 184 down to the stability line, exploring thereby the much studied 
"island" of stability. Here our barriers are in reasonable agreement with other 
microscopic calculations that have been made in this region: this constitutes 
an additional check on the overall validity of our calculations, and there 
is certainly no indication of a tendency for our calculational procedure to 
overestimate barrier heights. 

Nevertheless, we have not yet discussed triaxiality, which could be ex- 
pected to lower fission barriers to some extent. Several studies of this ques- 
tion have been made, the most extensive probably being that of Ref. (see 
that paper for other references), in which the barriers of 15 heavy and su- 
perheavy nuclei, some close to the stability line, others highly neutron-rich, 
were investigated within the framework of the ETFSI method and found, 
when triaxiality was taken into account, to be lowered by 0.6 MeV on av- 
erage, with a maximum shift of 1.3 MeV. It was concluded that this effect 
is probably negligible, given the overall discrepancy between calculation and 
experiment (see Section 1). Nevertheless, two counter-examples in both of 
which triaxiality is claimed to lower the barriers by the enormous amount of 
4 MeV are to be noted: 310 126 [[^J and 258 Fm []3"5| . The latter case, which was 



regrettably overlooked in Ref. Po| , is particularly disturbing, since it occurs 
in a region of the nuclear chart where the barriers of neighbouring nuclei are 
much less sensitive to triaxiality. The former of these two nuclei is too heavy 
for us to check with our codes, but we have repeated the latter case using 
the ETFSI method, and find that triaxiality lowers the barrier by 1.6 MeV. 



This is a significantly smaller effect than the 4.1 MeV claimed in Ref. [36|, 
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but it is still the largest triaxiality shift that we have found, and really too 
large to be neglected. One is thus forced to the conclusion that while in most 
cases triaxiality can reasonably be neglected, as far as barriers are concerned, 
there will be isolated cases where this is not possible. The situation is most 
unsatisfactory, since such cases can only be identified by first performing the 
full triaxial calculation, and doing this for the nearly 2000 nuclei considered 
here would be prohibitively time-consuming. It seems that the best that one 
can do at the present time is to proceed as we have done here, neglecting 
triaxiality, and bear in mind that a few of our barriers, including some low 
enough to be relevant to the r-process, will be considerably overestimated 
because of this approximation. However, the barriers of our proton-deficient 
nuclei in the vicinity of N — 184 will remain much higher than previously 
believed. 
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